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The two-particle problem within a nonequilibrium many-particle system is investigated in the 
framework of real-time Green's functions. Starting from the dynamically screened ladder approx- 
imation of the nonequlibrium Bethe-Salpeter equation, a nonequilibrium Dyson equation is given 
for two-time two-particle Green's functions. Thereby the well-known Kadanoff-Baym equations are 
iy-} ' generalized to the case of two-particle functions. The two-time structure of the equations is achieved 

£SJ ■ in an exact way avoiding the so-called Shindo approximation. For the case of thermodynamic equi- 

librium, the differences to former results obtained for the effective two-particle hamiltonian are 
discussed. 
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I. INTRODUCTION 
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This paper is devoted to the kinetic theory of many-particle systems which are able to form bound states. To 

• i-H i be specific, we will consider the case of two-particle bound states, i.e. the bound states are to be thought of, e.g., 

as hydrogen-like atoms or ions. Our focus will be the derivation of a kinetic equation for the distribution function 
of the (possibly excited) bound states. Starting with papers of, for instance, Waldmann Snider and Lowry 
, McLennan and Lagan ^|-|(| , Klimontovich and Kremp Q] , the derivation of kinetic equations has appealed great 
l — 1 j interest over many years. Usually (see, e.g. the review article of Klimontovich et al. H) the two-particle density matrix 
is split up in different parts with respect to some projection operator which projects onto the space of the bound 
states. Often the states are taken to be those of the isolated atom. The diagonal matrix elements are considered 
. to be the distribution function for the respective bound state. However, in a dense system, it is not clear if the 
diagonalization of the density matrix with respect to the unperturbed two-particle states is a good approximation. 

It is well known that many-particle effects like dynamical screening, self-energies or phase space occupation may 
have an influence on the two-particle properties. A unique description of these effects within the investigation of 
nonequilibrium behaviour can be given in the framework of the real-time Green's functions technique. For the single- 
particle functions, this can be done with the Kadanoff-Baym equations for the correlation functions ga ■ In this paper 
"--■^ , we aim on the derivation of similar equations on the two-particle level. 

Some remarks on the bound state problem in equilibrium seem to be useful first. The investigation of bound states 

• i-H | in dense systems, which are in the state of thermodynamic equilibrium, has been the topic of a lot of papers, see the 

references cited in the monographs || and PQ] . In the framework of the Green's functions method, a proper starting 
point is the so-called Bethe-Salpeter equation for the two-particle causal Green's function |ll| 



5ob (12,l'2')=Sa(l,l')5fc(2,2') (1) 



•l— I . r 

X ' +i did2dMg a {l,l)g b (2,2) K ab (12,i2) g ab {l%l'2') . 

. 

The kernel of this integral equation K ab is a four point function. In a diagrammatic language K ab comprises all 
irreducible diagrams in the particle-particle channel. The effective interaction K ab has a dynamical character. This 
makes the structure complicate: although one is interested only in the two-particle Green's function in the particle- 
particle channel with t\ — £2 and t[ = t' 2 , by the dynamical kernel K the knowledge of a Green's function with three 
times is enforced in the integral term. In Fourier space (or within the Matsubara technique) , this is adequate to the 
problem that, for the determination of the two-particle Green's function dependent on one frequency, a more general 
function dependent on two frequencies has to be known. A way out of this dilemma has been tried by applying 
the Shindo approximation jlij in which the causal quantity with two frequencies is constructed from that with one 
frequency. Then one gets a closed equation for the causal Green's function. There are few estimations on the range 
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of validity of this approximation. It is an exact relation for a static interaction K. It has been argued that the 
Shindo approximation reflects a first order approximation with respect to the retardation of the effective interaction 

In order to evaluate the Bethe-Salpeter equation, one has, of course, to use some approximation for the effective 
interaction K ab . In the simplest one, which has a dynamical character, K ab is just the dynamically screened potential 
V^f . We will come back to this point later. 

With help of this Bethe-Salpeter equation, some important questions concerning the properties of two-particle states 
in a plasma could be discussed. An effective Schrodinger equation was derived which has some important corrections 
in comparison with that for an isolated atom: (i) phase space occupation factors, (ii) exchange self-energies (Hartree- 
Fock), (iii) a dynamically screened effective potential and (iv) dynamical single-particle self-energy corrections. It has 
been shown that for localized states there is to a large extent a compensation between the effects (i) and (iii) on one 
side and (ii) and (iv) on the other. It follows that the binding energies of (at least the low lying) bound states are 
not changed considerably in comparison with the isolated atom. In contrast there is a large shift of the continuum 
edge by the self-energy corrections. This results in a lowering of ionization energies with increasing plasma density 
and leads, at the end, to the well-known Mott effect. An effective wave equation was solved numerically in |jT^,|l8[, 
for a discussion of the result see Kraeft [ p9| . 

However, the results obtained for the effective Hamiltonian have also some serious shortcomings. There occurs a 
division by Pauli-blocking factors 1 — f a — f b what causes spurious pole structures for highly degenerate systems. 
Further, the effective Hamiltonian has static contributions which lack a clear physical interpretation (see e.g. p(i|). 
These static parts vanish for a nondegenerate system. 

An other approach has been given by Schuck and co-workers |2f|j2^ ]. They postulate that Dyson equations exist 
for two-time causal and retarded Green's functions, respectively. Expressions for the self-energy operator (also called 
mass operator or effective Hamiltonian) are derived by comparison with the respective equations of motion. It remains 
also unclear in this approach what approximation (if any) is connected with the assumption that such two-time Dyson 
equations for the investigated functions and the inverse of those functions, respectively, do exist. Also in this approach, 
there occours the problem of the division by Pauli-blocking factors. 

Looking at the complicate structure of the Bethe-Salpeter equation, there arises the question if a general formulation 
in terms of a two-time Dyson equation is possible. The approach using the Shindo approximation is per se an 
approximative one. Schuck et al. do not consider the Bethe-Salpeter equation at all but postulate just this two-time 
structure. There are some strong arguments in favour of the possibility that exact equations exist for two-time Green's 
functions which contain solely two-time quantities (the inner structure of a self-energy operator could be complicate, 
however). In a chemical picture, for instance, in which bound states are interpreted as a new species, one would 
expect equations for correlation functions g^ b (t,t') similar to the Kadanoff-Baym equations in the single-particle case 
with some two-particle self-energy functions. The spectral information, e.g. binding energies, damping etc., should 
follow from the retarded two-particle Green's function g^ b (t,t'). The correlation function g^ b (t,t') is for t = t' just 
the two-particle density matrix. 

In nonequilibrium one can derive an equation of the same structure like in Eq. (|l]), however, the time integrations 
have to be performed then on the Keldysh contour. The functions involved in the integral term consist of many 
correlation functions because they are dependent on three times at least. Schafer et al. Q considered the dynamically 
screened ladder approximation for the polarizability in a semiconductor within the Keldysh formalism. They gave a 
formulation for functions depending on three times or - after Fourier transformation - on two frequencies. At the 
end, however, they used the Shindo approximation for these two-frequency quantities in order to get kinetic equations 
for single- frequency functions. 

There were attempts to generalize the Shindo approximation to functions in the time domain p4| . It was also 
tried, ]2{| , to generalize the approach of Schuck et al. within the nonequilibrium real-time Green's functions method 
postulating a Dyson equation for the retarded function gf lb {t,t') — 8(t — t')(g> b — g^ b )- In both approaches similar 
results were achieved. The equilibrium results could be reproduced. Thus the same shortcomings arise for degenerate 
systems. 

We will present here a new approach [p6| to this problem within the real-time Green's functions method. Thus 
we are able to describe nonequlibrium systems. Results for thermodynamic equlibium will appear as special case of 
the more general equations. In this first part the nonequilibrium Bethe-Salpeter equation is considered in a concrete 
approximation, the so-called dynamically screened ladder equation. This is the simplest approximation in which the 
effective interaction is of dynamical nature. This will enable us to identify the underlying algebraic structures and to 
keep the equations as simple as possible. The general scheme will be investigated in a subsequent paper |27[ . 

The structure of this paper is as follows. In Sec. |l| the scheme of the real-time Green's function method for single- 
particle Green's functions is summarized. The difficulties of the Bethe-Salpeter equation are discussed in Sec. III. 

The dynamically screened ladder approximation is considered in Sec. IV. The Bethe-Salpeter equation is written 
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down in this approximation for two-time functions on the Keldysh double time contour. The first terms of the 
pertubation expansion with respect to the dynamically screened potential are evaluated. The two-time structure will 
be achieved by applying the semi-group property of the ideal propagators. After that, the algebraic structures turn out 
to be similar to those of the non-equilibrium Dyson-Keldysh equation in the single-particle case and the two-particle 
self-energy functions can be identified by comparison. The thermodynamic equilibrium case is considered in Sec. 
The structure of the two-particle self-energy which can be understood as an effective Hamiltonian is discussed. The 
results are compared with the former ones [ fi^|]l5|]lO|]2(| . It will turn out that only in the nondegenerate case and in 
the static limit, one is led to the same results. 



II. SINGLE-PARTICLE QUANTITIES 



Let us summarize shortly the scheme of real-time Green's function technique in the single-particle case. The 
equations are given on a double-time contour, on the so-called Keldysh contour ]2q]29[ |. Objects of the algebra are 

matrices of causal and anticausal Green's functions, g a and g a , and the correlation functions ga defined below 

ga(rit,r[t') = G(t - ^(r^v'J) + G(t> - ^(r^f) (2) 

= 0(t - 04<*.(r 1 ,t)*I(ri,« , )> + 0(f - t)(±)4<*t ( r ;,f )* B ( rijt )) , 
i i 

S a M,r'i<') = 9(f - 0^(n*,riO + - *)<?a (rii,r'iO • 

Here a is the species index, Vfft and ^> a are creation and annihilation operators of second quantization with the 
commutation relations 



[* o (r,i),* 6 (r',t)] T = 
[*t( r ,t),*t (r / >t)]:F = 0< 



(3) 



The upper sign (meaning the commutator) holds for Bosons and the lower one (anticommutator) for Fermions. Spin 
is not written explicitely here. 

One can see that these elements are not all independent. It turns out that the equations get a more convenient 
structure if one introduces two other quantities, and g£, defined by 



ff * = e(t-07<*a(r ls t)*i(/ 1 ,*') : F*i(ri,f)*-(ri,t)) 



O(i-i') g>(t,t')-g<(t,t') 



-6 (f - t) - (*„ (n , t)& a (ri , T *J (r[ , t')* a (n , t) ) 



Q(t'-t) g<(t,t')-g>(t,t') 



(4) 



(5) 



It follows that 



9a = 9a 
9a = 9a 



9a 
9a 



9a 
9a 



9 a , 
9a- 



(G) 



The nonequlilibrium Dyson equation on the Keldysh contour reads 

9 a (^ 1 ')=9 a0 Ci-,l')+ I dldl5 a0 (l,l)S Q (l,l)3 a (l,l'), 



(7) 



with 1 = Tx,t\ etc., and g ^ being the ideal functions and E a the self-energy. The time integrations are performed on 



the Keldysh contour, see Fig. [J. The underlined quantities are causal ones for both times on the upper branch of the 
contour, g a — and anticausal ones for both times on the lower branch, ~g a = g~~. If the first time is on the upper 
branch and the second one on the lower, one gets g^ = gt - Fixing the first time on the lower and the second time on 
the upper branch gives g> = g~ + ■ Working on the Keldysh contour has the advantage that well-developed schemes 
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of functional derivatives and diagrammatic techniques known from equilibrium ]30| , pl[ | can easily be generalized to 
nonequilibrium situations, see, e.g., [f^ , [j2"| -f34"| . 

Eq. (0) is equivalent to a set of four equations which are, however, not all independent. Therefore, it is often more 
convenient to work with the following form of the nonequlibrium Dyson equation for the correlation function g< 



oo oo 

9a(t,t')=g< (t,t') + J dh J dt 2 [g< (t,h)^(h,t 2 )g^(t 2 ,t') 

to to 

+g* (t, h) E< (t x , t 2 ) g£(t 2 ,t') + g* (t, h) £*(ti, t 2 ) g< (t 2 , *') 



(8) 



which has to be supplemented by an equation for g a 



R/A 



CO oo 



9a(t,t')=g« (t,t')+ dh dt 2 g^ (hh)^(h,t 2 )g^(t 2 ,t'). 



0) 



TO CO 



Here only the functions's dependence on the times was written explicitely in order to save space. Often the initial 
time is considered in the limit t — > — oo. The quantity g^ is connected with g^ by hermitean conjugation 



g£{Tt,T>t>) = [g*{T't\vt)] 
An other form of Eq. (0) one can find is [p2[ 



oo oo 



Sa(M') = sf(Mo) <£(*o,to) 9?(to,t) + / dh / dt 2 gf(t,h)Z<(h,t 2 )g*(t 2 ,t) 



to to 



The corresponding differential equations read 



r d 


V 2 




2m c 


" d 


V 2 


l dt- 


2m ( 



E«(t,t) ff <(t,0 + E<(t,t)G^(F,f) 



(10) 
(11) 

(12) 
(13) 



These are the well-known generalized kinetic equations (Kadanoff-Baym equations) for the single-particle functions. 

In the following chapters the one-particle self-energy will be needed in a special approximation which is called 
^^-approximation. Here one has 

E a (l, 1') = Ef + *Zf (l, 1') g a (l, 1') , (14) 

with T, H being the Hartree self-energy. As for the Green's functions, cf. Eqs. ^^), there is a set of functions describing 
the dynamically screened interaction V_ ab 



vg{ht') = v ab s(t - o + e(t - t')vg>{t,t') + e(f - t)v a s b <(ht>) , 
v s ab (t,t') = - + e(t - t')vg<{t,t') + e(t' - t)v; s fc >(M') , 

V a s b R (t, t') = K 6 *(* - f ) + 9(* - H) [Vf b >(t, - 7/ 6 <(t, f)" , 

Kf/(i, = K 6 *(* - + e(f - t) kf b <(t, f ) - v;f>(t, f )' • 



(15) 



Here the correlation functions are defined by 



V* b <(t,t')=J2V ac Lf d (ht')V db , 



(16) 



where are the correlation functions of density fluctuations 

* L afc( r i t i' r 2*2) = (Spa{ri,h)Sp b (r 2l t 2 )) , 
^o6( r i*i' r 2*2) = (8pb(r 2 ,t 2 )5p a (ri,h)} , 



with <5p a (r,i) = *t(r,f)¥ a (r,i) - (*t ( r , t )$ ( r , t )). 
It follows that, e.g., V* <(1, T) = V£>(1', 1) and Vg(l, 1') 



) = 1), but V* «(1, 1') = V& A (1', 1) 



(17) 
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III. BETHE-SALPETER EQUATION 



The two-particle Green's function is determined by the so-called Bethe-Salpeter equation 

. 9ab (12,l'2')=.g a (l,l').g 6 (2,2') (18) 



+i J dld2dld2g a (l,l)g b (2,2)K ab (12,12)g ab (12,l'2' 



in which by introduction of the effective interaction kernel K ab formally a closed equation is achieved for the four-point 
function. Here, the kernel K ab is the sum of all diagrams irreducible with respect to a cutting of two single-particle 
lines. The Bethe-Salpeter equation can be understood to hold in various contexts: for T = 0, for the imaginary-time 
equilibrium Green's functions in the Matsubara technique, or for the real-time Green's functions on the Keldysh 
contour. 

The properties of a pair of particles should follow from this equation in the so-called particle-particle channel. If 
one considers the causal two-particle Green's function in this channel (t\ =t%=t] t[ = t' 2 = £'), one has 

i 2 ga b (rir 2 t,r' ir ' 2 t') = 0(i - t')(*„(r 1 ,t)* 6 (r 2) t)*J(r' 2 ,0*i(ri,0) (19) 
+ 6(t> f )(*t (r[, f , 0* 6 (r a , main, t)\ 

= e(t-ty 9 > b + 9(t'-t^ 9 <. 



On the right hand side of Eq. (|18|), however, there occurs a function depending on three times t%, t 2 , and t' what is 
enforced by the dynamical character of K ab . This Green's function consists of six different correlation functions 

(20) 



= 0(tl 


- h)0(i 2 


-t')(* (f lj t 1 )* 6 (f a ,i a )*t(^, 0*3,(^.0 


± 0(h 


-h)0(h 


- t')(M^, f 2 )*«(f i, hMfo, 0*H(ri, 


+ 8(h 


- t')9(t' - 


-f 2 ) ( * Q (f 1 , t ! ) tf b (r 2 , t' ) *t (ri , f ) (r 2 , t 2 ) | 


+ 9{i 2 


- t')9(t' - 


-fx)(*6(r 2 ,f2)* t o (r' 1 ,t , )^(r 2 ,0*a(ri,fi)] 


± 6(t' 


- ti)0(ti 


- f 2 ) (ri , t')^(r 2 , t')*a(ri , fi ) * b (r 2 , f 2 ) 


+ 6(t' 


~ h)0(h 


- f x ) ( *t (ri , (r 2 , t')*6 (f 2 , f 2 )¥ a (f ! , fx ) 



Only a static interaction in ([l8]) would enforce fx = ^2, and the function would turn into the two-time causal one jis|). 

In principle, one could of course try to solve the Bethe-Salpeter equation for a function depending on three times 
(equivalent to a function depending on two frequencies) and then extract from this the information one is interested 
in. This, however seems to be to complicated. In a number of papers it was tried therefore to work with equations 
which involve exclusively two-time functions. This was achieved in two ways. The first approach ]l^ , [r^ | uses within 
the Matsubara technique the so-called Shindo approximation in which the two-frequency function is constructed from 
the single-frequency causal Green's function. This is possible in an exact way for a statical interaction and therefore 
it was argued that for arbitrary K ab this would be correct in first order with respect to the retardation. In the other 
approach plj ], a closed equation for the causal two-time or single-freqency, respectively, function is postulated. The 
effective hamiltonian (two-particle self-energy) is determined then by comparison with equations of motion. 

In any case there are closed equations for combinations of the correlation functions g< b defined in ( |l9| ) . But this 
means that the information contained in the correlation functions a^l^^ b) and (^b^a^l^a) (third and fourth 
term in Eq. is neglected. Therefore such a closed equation for the causal two-time Green's function can exist only 
in an approximate way. 

In the next section we will show that closed equations exist for other combinations of correlation functions. 

IV. DYNAMICALLY SCREENED LADDER APPROXIMATION 

In analogy to the single-particle case, one is interested to get information on the statistical properties, carried by 
the two-particle density matrix, as well as information on the two-particle dynamics, the spectral information. Below 
we will see that it is not a trivial question to say which quantity carries this information. 
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One of the quantities of interest is the following two-time correlation function 



g< b {v 1 v 2 t,v' 1 v' 2 t') = i(*t (r' 1; t')^ b (r' 2 , t')* b (r 2 , t)9 a (r l: t) 



(21) 



In the case t' = t, the quantity i 2 g< b is just the two-particle density matrix /^(rii^r'^, t ). 

We will use in the following the real-time Green's function technique in the Keldysh formulation. In the present 
section the Bethe-Salpeter equation will be considered in a concrete approximation: for the effective interaction kernel 
K a b, the dynamically screened potential V^ b is taken. Then one has on the Keldysh contour 



9^ b (ht 2 ,t' 1 t' 2 ) = g a (t 1 ,t[)g i (t 2> t' 2 ) + I dhdhgJt 1 J 1 )g b (t 2 ,t 2 )iVZ b (hj2)gJht 2 ,t / 1 t' 2 ) . (22) 



Iteration of this integral equation leads to ladder-type terms. 

This Bethe-Salpetcr equation ( p2] ) will be considered in the following in the special case t\=t 2 =t and t[ =t' 2 = t' . 

A. Two-particle Green's functions depending on two times 

Although the physical times t\ and t 2 as well as t[ and t 2 are equal, there are still 16 possibilities to fix the times 
ti, t 2 , t[, and t' 2 on the upper and lower branches of the Keldysh contour. Fixing the times t\ and t 2 on the upper and 
t'i and t' 2 on the lower branch of the Keldysh contour, one gets an equation for the correlation function g^ b = g^ b 
defined in Eq. (pl|). Below there are given three other important cases. 



9ab 



9ab 



9ab 



-(r 1 r 2 t,r' 1 r^') = e(t-i')-j(*i(r' 1 ,t')* b (r 2 ,i)* Q (r 1 ,t)vI/t(r' 2 ,0 

+0(f - t)i(* 6 (r 2 ,t)*t (ri.iO*^, tf)*a(r ls t)) , 
( ri r 2 t, v\v' 2 t' ) = 6(i - 04(*6 (r' 2 , t')^a(ri, *)* 5 (r 3 , (r[, t') 



+0(f - t) (*„ (n , m\ (r 2 , f )*J (ri , f , t) 



1 



( ri r 2 *,rir 2 = * B (n, *)* 6 (r 2 , i)^(r 2 , 0*tW, ) • 



(23) 



(24) 



(25) 



To give an example, the time-ordering on the Keldysh contour is shown for g^ b ' + in Fig. [|. All other functions 

g" b nS with the greek indices equal to "+" or "— " can be expressed in terms of the six correlation functions involved 
in Eqs. © and (§§|2|). 

We define the following retarded and advanced quantities 



e{t-t)i(g+ + <- -g+->-+-g^+- 
e(t-f)j([*i('i.*'), [^(r 2 ,0.*6(r 2! t)*a(ri,t) 



(26) 



G^r^i, r[r' 2 t' ) = 0(f - t)(-t) - 5 + 



Sab 



Safe 



(27) 



e(t'-t)—( 



*a(ri,t), 



* fe (r 2 ,t),^(r 2 ,t')^(ri,t') 



In order to achieve the nested commutator structures, it was used that operators with equal times can be interchanged 
according to (|^). Interestingly enough these nested structures were found also by Rajagopal and Majumdar in their 
analysis of double dispersion relations for the two- frequency causal Matsubara Green's function (Appendix II of 
We show in Appendix ^ as an example how the function G^ b is connected in thermodynamic equilibrium with the 
analytic continuation of the two-frequency Matsubara Green's function. 

The functions G^ b A have the following properties: 
(i) They are connected by hermitean conjugation 



G&farai.r^f) = G£ 6 (r , 1 r , a * , ,r 1 r a t) 



(28) 
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(ii) Both functions have the property of crossing symmetry, i.e. 



(29) 



(iii) The inhomogeneity in the equations of motion for these functions consists of 5-functions only (without Pauli- 
blocking terms). This is easy to see from Eqs. ( p6|]27| ). Derivation of the Heaviside-function gives a 5-function 6(t — t'), 
and the commutation relations for the field operators lead to [<5(ri — r' 1 )i5(r2 — r' 2 ) ± 5 ab 8(r x — r' 9 )8(r 2 — r'i)]- 

(iv) For vanishing interaction between the particles a and b, the correlation functions in Eqs. (^6 27) can be contracted 
into products of single-particle correlation functions g<, and one gets 



G&(rir 2 t, v' x r' 2 t' ) = i gf(r x t, r' x t')g?(r 2 t, r' 2 t') ± S ab i g£{r x t, v' 2 t')gf{v 2 t, v' x t') , 
with the retarded single-particle function g^ defined by Eq. (^). 

(v) The difference of the retarded and the advanced functions defines a spectral function 

A ab {nr 2 t, r' x r 2 t' ) = iG* b [r X T 2 t, v' x v' 2 t' ) - iG^nrrf, r[r' 2 t' ) 



(30) 



(31) 



9ab 



9ab 



9ab 



which gives in the case of equal times (t = t') 

A ab (v x r 2 t, r' x r' 2 t ) - 8(r x - r' x )S(r 2 - r' 2 ) ± 8 ab S(r x - r' 2 )5(r 2 - ri) . 



(32) 



B. Transformation of the Bethe-Salpeter equation into a Dyson equation 



In this subsection the dynamically screened ladder equation as a special approximation of the BSE is transformed 
into a Dyson equation in which the occuring two-particle Green's functions and two-particle self-energy functions are 



e _p< 

considered in the diagrammatic form shown in Fig. |2| It is analyzed first for g^ b ; details are presented for different 



dependent on two times only. For this purpose, the_perturbation expansion of the Bethe-Salpeter equation (f22|) is 
considered in the diagramma 
orders of Vf b in Appendix |A|. 



We search for (and, indeed, find) the following structures [cf. the corresponding equations for the single-particle 
functions, Eqs. || and || 



g< b = g< b + g* 

^ab = Sab + &a 



V ab 
V ab 



9ab + GabVab G ab + $a 



Van 



G 



ab 



G 



ab 



(33) 
(34) 



with the definitions Q< b (t,t') = g< (t, t') g< (t,f) , Q* b {t,t') = ig* Q (*,*') g* (t,f) , and # 6 (t,f) = 
(— i)gao{t, *') <?{To(*' *') ■ -^-11 quantities in the above equations depend on two times only. The integration of in- 
termediate times runs in the interval [to, oo] like in Eqs. (^|^). 

The zeroth, first and second orders for g< b with respect to the two-particle self-energy are given by 



<(0) 

9ab 

<(i) 

9ab 



g < 



b I 



Vab 



Sab + Gah^ahGab + Gf> 



ab ab^ ab 



* ab 



Vab 



'ab ) 



(35) 
(36) 



<(2) 
9ab 



Vab 



^ab 



Vab 



-'ab 



Gab + Si 



Vab 



-'ab 



rR _< r A 

yab a ab y a, 



Vab 



ab 



Vab + ^ b 



yah 



Vab 
Vab 



Gab + Sab a abGab 



Vab 



y a b 



y n h 



(37) 



The self-energy functions a^ b and E^ fe , respectively, are identified by comparison with the expansion terms of the 
ladder equation, Fig. g, then. 

All functions in the above equations ([53|) and (|34|) are understood to depend on two times. The key idea in order 
to achieve such a two-time structure of the equations is to use the semi-group properties of the ideal single-particle 
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propagators g^ and g£ (the time-local Hartree-Fock self-energy could also be included). In particular one has for 
any time i with t > i > t' the following relation 



(38) 



There is no time integration in the above equation. Analogeously, for the advanced function with t < t < t' one has 
(integration with respect to r 2 suppressed) 



For the ideal one-particle correlation functions g< , there follows 



(39) 



(40) 



9%(t,t') = igK (t,t)g%(i,t') for t > i, 
9%(t, if) = H) 9%(t, t) g£ Q (i, t') for i<t', 

9%(t,t') = (~i 2 )gK (tJ)g%(t,t)g£ (i,t') for t>t and t < if . 

Proceeding in the manner presented in Appendix ^ and comparing the results with the anticipated structure 
[Eqs. (^||),(|3j)], we get the following expression for the retarded self-energy function 



E* (f, = i£?(t, t') g« (t, t>) + i£*(t, ff « (t, f ) 

+*<$,(*, f) iVf b >(t, if) g* (t, t') + ig< Q (t, t') iV a s b R (t, f ) g* (t, t') 
+ig* (t, t') iV a s b <(t', t) g* (t, t') + ig< (t, t') iVf b A (t', t) g R (t, t') , 



(41) 



where the one-particle self-energies have to be used in first order of the dynamically screened potential V s [cf. 
Eq. ©],i.e. 



£*(f , f ) = E?(t)*(f - t') + ^>(t, g R (t, + 5 < (t, . 

The correlation function crS is found to be 

a<(t,t') = X<(t,t')g< (t,t') + X<(t,t')g< (t,t>) 

where the single-particle self-energy function £ a is given by 

£<(t ) t , )=*t£ < (t 1 t , )ff<(t,t / )- 

The diagrammatic structure of the two-particle self-energy functions is shown in Fig. ^. Primarily, these functions 
consist of naked lines because we worked in first order with respect to the dynamically screened potential. Howe ver, 
all diagrams necessary to dress the lines could be found in higher orders of the expansion, see Appendix A 2 and 
especially Fig. ||. 

The self-energy functions E^/" 4 and a^ b are functionals of single-particle Green's functions. That is the reason why 
it is sufficient to consider the two equations ([33]) and ( |34| ) in order to determine the correlation function g^ b . In higher 
approximations, the two-particle self-energy is expected to be a functional of two-particle correlation functions, too. 
Then one would need also the equations for the other three quantities defined in Eqs. ( p3f|25| ). The full scheme of 
equations reads 



(42) 
(43) 

(44) 



9ab — y a b t y a b 



Vab 



-"ah 



9ab + y a b a ab G ab + V ab 



Vab 



G 



ab 



with $ = {++, }; {+-, -+}; {-+, +-}; { , ++} , 



b t y a b 



G ab = <d t 
Gab = Sab + Sa 

For the see next subsection. 



Vab 



Vab 



g: 



i ■ 



u ab 



G 



ab " 



(45) 

(46) 
(47) 
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These equations are not all independent. G^ b and G^ b are connected by hermitean conjugation. Further they are 
linear combinations of the preceding four functions g f b ac cording to Eq. (|26j). The system of equations is consistent, 
i.e., combining the equations for the g% according to ( j26| ) one gets the Dyson equation d46|). 

Often it is more useful to consider differential equations 



i--H° 



Vat 



(It 



i) g* b (t, t') + ff * (t, t) G*(i, t<) 



(48) 



There are additional equations for the propagator functions G^ b and GA. As these two functions are connected by 
Hermitean conjugation, only the respective equation for G^ b is written down. 



dt a 



Vol 



G B ab {t : = 5(t-t')+ diX*(t, t) G* b (t, f) 



(49) 



The Eqs. @-^) and (||-||), respectively, can be considered as the most important result of the present paper. 
The latter equations are the two-particle counterpart to the Kadanoff-Baym equations in the single-particle case. 
Thus, these equations are the proper basis for the description of two-particle properties. 



C. Algebraic structure of the two-particle self-energy 



Similarily to the four different two-particle Green's functions considered in Eqs. ( |2l| , |23j-|25j) self-energy functions 
<T ah' 5 ' 7 ,l5 ( r i r 2^ r i r 2^' ) w ere introduced in (f45|). Within the present approximation they are given explicitely below . 
The function cr+ h + -~ ist just given by a< b of Eq. (||). Making use of V^l, T) = V^<(1', 1), we get 



'ab 



,- + ■ + - 



V^<(nt,rif) + V£ <(r 2 ty 2 f) + V a s b <{v lt y 2 t>) + V^<(r 2 *,riO 
i [V a s a <(r lt , rtf) + V b s b >(r 2 t, r 2 i') + V a s b (nt, v' 2 t') + Ff a (r 2 t, r[t') 

Xfl'o,o( r l*> r/ l*05^o( r 2*> r 2* / ) i 

* [\^ > (r 1 i, riO + Vjf <(r 2 t, r' 2 t') + Ff 6 (nt, r 2 i') + V^(r 2 i, r^')" 
Xffa,o(ri*>r' 1 t')5'6 < :o( r 2*) r 2*')! 



-,-H 



' a6 



The retarded two-particle self-energy is given then by 



(50) 
(51) 
(52) 
(53) 

(54) 



+Q(t-t')i 



'ab 



'ab 



' ab 



' ab 



The term which is local in time consists of the single- particle Hartree and Hartree-Fock self-energies as well as the 
Pauli-blocking contribution 



E^nr^r^t') = (nritWra - r' 2 ) + Ef* (r 2 r 2 t)«5( ri - r[) 



(55) 



+ 



ig<{v 1 t,T l 1 t)5{T 2 ^v' 2 ) + ig<(v 2 ty 2 t)5(r 2 -v l 2 ) V ab {r[ - r' 2 ) 



Inserting the expressions for the cr* b , Eqs. (50-53(), into Eq. |4|, one gets indeed Eq. pi] ) 
The advanced quantity is hermitean conjugated, and one has 

E^( ri r 2 t, r[r> 2 t' ) = E° 6 +( ri r 2 , r[r 2 t )5(t - t>) 
+Q{t'-t)(-i) 



(56) 



,T++> 

a ab ~ a ab 



, + + 



' ab 



with 
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(57) 



E°t( rir2 t, r'^t' ) = (nrit)*(r a - r' 2 ) + E^ (r 2 r 2 t)£( ri - r' x ) 

-r- 2 )[i fl <(nt,rit)5(r 2 - r 2 ) + ig<(r 2 t,r' 2 t)6(r 2 - r' 2 ) 

Often it is useful to consider Fourier transforms with respect to the difference time t — t' according to 

F(u, T) = J d(t - t'y^-^Fit - t', T) , (58) 

with the "macroscopic" time T = (t + 1')/2. In thermodynamic equilibrium there is no dependence on this time. 

The analytic properties of the two-particle self-energy are more involved than those of the single-particle self-energies 
|?5f because there can occur off-diagonal matrix elements already in the spatially homogeneous case. The two-particle 
self-energy E^ b (w,T) can be split into hermitean and antihermitean parts according to 



E3,(121'2>T) = E«(121'2>T) + E&(121'2',wT) 



(59) 



with E^(121'2',wT) = E^(l'2'12,u;T)* and E^,(121'2', ujT) = -E$,(l'2'12, ujT)*. These functions can be con- 
structed in the usual way. The are connected by Kramers-Kronig relations, see Appendix 0. 



V. THERMODYNAMIC EQUILIBRIUM 

A. Two-particle Dyson equation 

In thermodynamic equilibrium only the spectral properties have to be determined, i.e., only Eq. ( |49| ) for the two- 
particle propagator has to be considered. Its Fourier transform is 

n H° ab ~ V ab E° b - E« b c °" (fi)] <?*(£)) = 1 . (60) 

This equation can be called two-particle Dyson equation or, likewise, Bethe-Salpeter equation. We want to draw the 
reader's attention on the fact that this equation is given here for the function G^ b , whereas in earlier attempts it was 
tried to formulate such an equation for the causal two-particle GF, |T^j2l[], or a function g^ b = Q(t — t')(g > b — g^ b ), 
©• 

The static part of the two-particle self-energy in (60) is given by E° & = E^ F + E^ F + (N ab — l)V ab . The self-energy 
consists of two different types of terms: some do not contain an interaction between the particles a and b whereas the 
others do. The first terms are due to single-particle self-energies. 

For the correlation part of the two-particle self-energy, there holds the same distinction: E^ b corr (f2) consists of two 
contributions according to 

Efr T (^)=Af b (0) + ^f R (0), (61) 

where the first one is due to the one-particle self-energies, whereas the second one describes an effective interaction 
between particles a and b in the many-particle system. We have 

oo 

A5(piP2,PiP 2> n)= / — -—3 — / 7^3^ — -= 62 

J (2^) J ( 27r ) " - w l — ^2 - ^3 + «0 

— 00 

x{[*Ka(<l> w l)*0a (Pi +q>^2) -iV;<(q,wi)^<(pi +q,w 2 )][«.g b > (p2,^3) -«ff b < (P2,^3)] 
+ [^(q,wi)i5 fc > (p 2 -q,w 2 ) -iV^(q,wi)i ff <(p 2 - q, w 2 )][ig>(pi, w 3 ) - ig< (p u 
x(2n) 6 S(p 1 ~p' 1 )S(p 2 -p / 2 ) 7 



and 



Kf R (PiP2,PiP 2 ,^) (63) 

OO 



(2^) 3 ,5(p 1 +p 2 -p' 1 -p 2 ) f 



(2ir) 3 £1 — uj\ — lu 2 — CJ3 -|- iO 
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X j{[«5a (Pl.^l) ~Wa (Pl,Wi)]iV^(pi -pi,-W 2 ) 

+*3a (Pi,wi)] [iV^(pi -p'n-wa) -^(Pi -P'i,-W2)]} [*sf (P^) -^(P^^s 
+ {[*5 6 > (P2 ! wi) -ig b < (p 2 ,wi)]iV a ^(pi -pi,w 2 ) 

+*5 fc < (P2,^i)] [iV^(pi -pi,u; 2 ) -iV t(pi -Pi,w 2 )]} [i5>(pi,w 3 ) -«5a (Pi> w s)] 

This can be further evaluated using a quasiparticle approximation 

±«5a (P>^) = 2tt(5(cj - e a (p))f a {uj) 
ig>(p,uj) = 2tt5(u - e a (p)) [1 ± /„(w)] , 



(64) 



with f a {uS) — {exp[/3(w — fi a ] =p 1] 1 . The correlation functions c an be expressed in terms of the dielectric function 
e(uj) and Bose functions n B (uj) — [exp(/3cj) — according to ]36| . |l0[ 



iV<(q,uj) = -2V ab (q)lme R (q )W )n B (w), 
iV^(q,w) = -2K b (g)Im £ ii " 1 (q ! w) [l + n B (w)] 

For the function A^ 6 (f2) there follows 



(65) 



A« b ( Pl p 2 ,p' lP2 ,0) = (27r) 6 <5 Pl , p; <5 p2 , p 







1 (2^) 3 J 


7T 



Im£ fi (q,wi) 



(66) 



1 ± / a (pi +q) +n B (wi) 



(27r) 6 <5 piiPi (5 P2:Pi 



- wi - e a (pi + q) - £b(P2) +«0 

a5,( p ' iP2 ,o). 



+ (a 



6, 1 



In the last line of this equation, the function A^ ? h (p' 1 p 2 , SI) was introduced which is the sum of the single-particle 
self-energies (in V s approximation) to be taken off-shell 



A« (p' lP2) SI) = E« pi, n - e b (p' 2 ) + Ef p 2! n - e a (pi 



(67) 



The other contribution to the two-particle self-energy, V^f , is given by 

oo 

V a f R (p 1 p 2 ,p' lP2 ,fi) = (2^) 3 J pi+P2 ^ i+p , | ^ [-Ime^Cpx-pi,^) 

— OO 

i±/.(pO + »*(«!) + (a ^ 6)1 ^ 2) 



x|k6(pi - pi! 



n - wi - e (pi) - £b(p 2 ) +*o 



(68) 



These two contributions to the two-particle self-energy look very similar. Replacing the Coulomb potential V ab by 
z a z b V with z a , z b being the charge numbers, one can see that for particles attracting each other, there is a compensation 
between these two functions. This is especially to be seen considering the functions integrated with respect to pi and 

P2 



d 3 pi f d 3 p 



(2tt) 3 J (2tt) 3 

OO 

/ \ f d 3 q f dcji 

y Za+Zb )J 1 — 



A^( Pl p 2 ,pip 2 , SI) + V a f K ( PlP2 , p' lP2 , SI) 



-lme R (q,wi) 



(69) 



1 fi-cJi -e a (pi +q) -e 6 (p 2 )+z0 



6, 1 
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In the case of a symmetrical plasma, z a — —z b , the right hand side of the above equation vanishes and it follows 

A« 6 (piPa.n) = -/ |^ / ^^f R (PiP2,PiP^^). (70) 

These expression for the two-particle self-energy have to be compared with the results of former papers fll2] , |l5| |. 
The notations are slightly different in comparison with ours, so one should compare the expressions of the effective 
hamiltonians. The total hamiltonian is H® b + V ab + H^Jui) with H® b + V ab being the hamiltonian of the isolated pair 

of particles whereas the medium-dependent part of the Hamiltonian is denoted by H^ b {uj). In the present paper this 
latter quantity is given by 

fl^HESa + ESTM (71) 
= E« F + E« F + N ab V ab - V ab + A R b (Sl) + Vg R (fi) , 

where w is to be understood as an parameter. 

The differences consist in the following: (i) now there a no additional static parts beyond the Hartree-Fock level, 
(ii) no division by Pauli-blocking terms occurs. Both things seem to be produced artificially by adopting a closed 
equation for the wrong quantity. 

On the other hand, the same results as before are received in the nondegenerate case. 



B. Limiting cases 



It is interesting to study some limitin g ca ses of our expressions. First, in the nondegenerate case the one-particle 
distribution functions in Eqs. (|66| ) and ( p8| ) can be neglected. The result is in aggreement with the nondegenerate 
limit of the former approaches using the Shindo approximation J^,[TJ|,|lo|,^o| . 

The second important limiting case is that of statical screening. Following Zimmermann |20[| we consider the case 
that the excitation energy into a pair of two free particles, e a (pi) + £b{p' 2 ) ~~ SI, is small in comparison with the energy 
lo\ occuring in the dielectric function. This could be a reasonable approximation for excited states. Then we have 



A? 6 (p lPa ,pil4 SI) = (2tt) 6 6(pi ~ p'J <5( P2 - p' 2 ) 



x^ aa (q) 1±/a(Pl+q) : nBM +(a 



-u>\ +i0 



(2TT) e 6( Pl -p' 1 )S(p 2 -p' 2 ) 



d 3 q 
(2tt) 3 

6, 1 «— , 



doji ' 

7T 



lme R (q, u\) 



2) 



d 3 q 
(2^ 



V a S Jq, Sl = 0)- V aa (q) ± / Q ( Pl + q) 



V£ (q, n = 0) - V bb (q)} \ ± f b (p 2 + q) + J 



(72) 



and for the effective interaction term 

V a f R ( PlP2 , p[p' 2 , SI) = (2ir) 3 S( Pl +p 2 - pi - p' 2 



duj\ 

IT 



lme R (pi-p'i,wi) 



T/ , M l±/a(Pl)+"£g(^l) , i 
V ab (pi -pj — + (a 



b, 1 



2) 



= (2tt) 3 s( Pl + P2 - pi - P2 ) v a s b ( Pl - pi, n = o) - v ab ( Pl - pi) i ± / a ( Pl ) ± h( P2 ) 



(73) 



Here it was used that Ime 1 is an odd function and that the even part of the Bose function rig (a;) is —1/2. Further 
one has 



vg*{q,n)=v ab {q)\i- j 6 ^ 



du>i Ime 1 (q, Wi ) 
SI — lui + iO 



(74) 



The terms in Eq. ( |72[ ) containing single-particle distribution functions and Eq. ( |73| ) can be combined with 
and N ab V ab , respectively, to give functionals and N ab Vf b (Sl = 0) of the screened potentials in the 

static limit. The remaining terms in Eq. (|7^) give a constant term. 
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The Dyson equation, Eq. (60), can be written then in the following form 

n - ] g* (n) = 1 , 

with the effective plasma Hamilton operator 



H°f = H° b +J2x? F \v s (n = o)\ 

c=a,b ^ > 



+ N ab v a s b (n = o) 



(75) 



(76) 



5 E 



d 3 r 



2 ^ y (2tt) 3 

c=a,b J y ' 



V c S c {q,Q)-V cc {q) 



Considering this effective Hamiltonian for a nondegenerate system, one can write (with V ab — > z a z b V) 



H C J = H a ab + z a z b V{r) + z a z b \V S (r, O = 0)- f (r) 

1 r " ' [y s (o,n = o)-y(o) 



(77) 



2 i 2 



For a symmetrical plasma with z a = — z b = 1 this leads to 

= H° ab - V(r) - ( [V s (r, Q = 0) - F(r)l - \v s (0, Q = 0) - V(0) 



(78) 



i.e., at small distances the last two terms compensate each other to a large extent and the interaction is given by the 
Coulomb potential 



hm Lff Q t = H° ab - V(r) + 0(r 2 ) . 

r— >u 

Adopting for V s (r, = 0) the statically screened Debye potential V D (r) , 

- 2 47re 2 

g 2 + k 2 



V D (r) = — exp (—rer) ; U%) 



(79) 



(80) 



one gets for the Hamiltonian in ( |77| ) 

flSf = H "b + *aZ b V D (r) - l - (4 + zf) ne 2 . (81) 

The two last terms combined give the well-known effective potential of Ecker-Weizel type |37j which has been used 
frequently in order to determine energies and wave functions of bound states in a plasma environment |l5|,|l^,[58],[39| . 



C. Effective wave equation and two- particle energies 

The Eq. ([30]) which determines the two-particle propagator G^ b was written down in an operator form. Using a 
representation one gets a matrix equation. In order to solve this equation it is favorable to use a representation in 
which diagonal elements are the main contribution. We follow here Kilimann et al. [ fl5|| , however, now it is not possible 

±1/2 

to achieve symmetric real and imaginary parts of the effective hamiltonian simply by multiplying with factors N ab 
Therefore the hamiltonian is split into hermitean and antihermitean parts, see Appendix H. This leads to 



(82) 



with H^Q) = H a + H b + V a b + E"^ + S^ rrW (0). The eigenvalue problem of the hermitean part of the hamiltonian 
reads 



H^ b (rt) nP,n) = E nP (n) n p,n 



(83) 
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The eigenstates nP, SI) where nP denote the quantum numbers and u> is a real parameter can be used as a orthonormal 

basis. The eigenvalues E n p(Sl) of this effective Schr"odinger equation are not yet the spectrum of two-particle excitions 
p5| . The latter follows from the spectral function A a b. 

In the representation with respect to the eigenstates \nP, Sl^, Eq. ^ reads (conservation of center-of-mass momen- 
tum already taken into account) 

[n - E nP (Sl)] G R nr , ft, 0) - Yl H nm(P, nA> ft = <W ■ (84) 

m 

In the following it is assumed that nondiagonal matrix elements of the antihermitcan part of the effective hamiltonian 
are small. Then the Eq. ^ has the approximate solution 

G^n' ft ft) = Q + i o-E n p{n)+iT nn {PQ) (85) 

-iT nn ,(PSl)(l-6 nn ,) 



[n + io - E nP (si) + iT nn (pn)][n + io- E n , P (n) + ir n , n ,(psi)] ' 

where it was introduced T nn / (P, SI) = iH^ n , (P, SY). For the coherent part of the spectral function there follows 

A nn {P, SI) = — — / M2 , ^2 1 nn\ • (^6) 

[si - E nP (si)Y + r£„(Pfi) 

According to this equation, the spectrum of the two-particle excitations is given by the roots E n p of 

Q = E nP (Sl) (87) 
whereas the damping is given by T nn (P,E n p) Jig ], 



VI. SUMMARY AND CONCLUSION 



Starting from the nonequilibrium Bethe-Salpeter equation in the dynamically screened approximation, we have de- 
rived a set of nonequilibrium Dyson equations for two-time two-particle correlation functions. The two-time structure 
of these equations was achieved in an exact way using the semi-group properties of the ideal one-particle Green's 
functions. The prize one has to pay for this simpler structure of the equation is that the two-particle self-energy in the 
Dyson equation consists now of irreducible diagrams in all orders with respect to the dynamically screened potential 
(in some sence this is similar to the transition from Feynman diagrams to Goldstone diagrams Irreducibility 
means here that a diagram cannot be cut with respect to a pair of single-particle lines which begin at equal times and 
end at equal times, i.e., two ore more interaction potentials have some overlap in time. 

For the further considerations we have restricted ourselves to a two-particle self-energy in first order with respect to 
the screened potential. The algebraic structure of the equations is not affacted by this approximation. It was shown 
that there is a set of equations for four two-time correlation functions. This generalizes the couple of Kadanoff-Baym 

equations for the one-particle correlation functions g^ (<?_| and g |_, respectively). In analogy to the single-particle 

case there is no closed equation for the correlation functions but always a coupling to other correlation functions. Only 
for two certain functions G^ A there exist closed equations. Thus, these functions are the two-particle generalization 

of the retarded (advanced) commutator Green's functions g^ A in the single-particle case and just these functions 
describe the propagation of a pair of particles in the nonequilibrium many-particle system. 

The case of thermodynamic equilibrium was considered in some detail in order to show the differences to former 
approaches. In former attempts |4(]|lJ] there were anticipated closed equations for the causal two-time two-particle 
Green's function. These equations were enforced by the Shindo approximation. The expressions for the effective 
hamiltonian are the same as those of the present paper only for the case of a nondegenerate system. The agreement in 
this special case is easy to understand taking into account that the difference between the used functions is of higher 
order in the density. 

For arbitrary degeneracy there are clear differences between the former results and ours. In the present results 
there is no division by Pauli-blocking terms. The only intrinsic static contributions of the effective Hamiltonian (the 
two-particle self-energy) are the Hartree-Fock single-particle self-energies and the Pauli-blocked basic potential. 

We can conclude that the proper generalization of the Kadanoff-Baym equations for two-particle functions is given 
by the system of equations ( p8| ) and (f49|). The algebraic structure of these equations was identified starting from a 
concrete approximation, the dynamically screened ladder equation. More general considerations how the self-energy 
functions can be determined in higher approximations will be presented in a subsequent paper p?j . 
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APPENDIX A: EVALUATION OF DYNAMICALLY SCREENED LADDER TERMS 

The aim of this appendix is to show the evaluation of the lowest-order terms in the dynamically screened ladder 
equation. Single-particle self-energy contributions and interaction terms have to be treated on equal footing. Special 
attention is paid to the transformation into a structure involving two-particle quantities which depend on two times 

Similar considerations are possible 

and 9ab 



only. The analysis is made here for the expansion of the function g< b = g^ b + 



for the other three functions gj b ' , g ab 



-.+ — + „-++- 



This is scetched in subsection A 3 



1. First-order contributions 



There are three diagrams of first order with respect to the dynamically screened interaction Vf b , see Fig. ^. Two 
terms have single-particle self-energy insertions of particles a and b, respectively. The third one is a ladder diagram 
with one rung. 

The first term with a self-energy insertion for particle a is given simply by [cf. Eq. 0] 



l[ 1} (t,t') = / dtxdh [ff<o(*.*i)^(tiJi)^o(*i.O+55o(*.*i)S<(t 1 ,ti)^ (*i,*') 



+9a,o(t, h) X*(ti,h) g< (h,t') 9b,o(t> • 



(Al) 



to 



In order to achieve the anticipated structure, one can use the semi-group properties for the correlation function 
g bQ (t,t!). In the first term on the right hand side, for instance, there holds t\ < i\ < t' , enforced by the advanced 

functions i x ) g£ fl ($u *')• For this case we can use in Eq. ( |AT| ) g bfi {t,t') = S^ (t,*i) (-i)s£o(*i>*i) 

Treating the other two terms in a similar way one gets 



A (1) (M') = J dt 1 di 1 {g^ (t,t 1 )g< (t,t 1 ) [(-i)^(ti,fi)^ (*i.*i)] HKo(M') ff&fr.O 

to 

+*s£o(*. h) g£o(t, h) [s<(ti, h) gfoih, h)] H)s£o(*i>* / ) 0m>(*i> 

+ig*; (t,t 1 )g b B ; (t,t 1 ) [t£? ; (ii,ti)s£ (*i>*i)] ^ofr^'Ko (M')} ■ 

This fits into the structure 



(A2) 



&ab ^ab Gab + Gab °ab &ab + ^ab^ab^ab 



(A3) 



(2) 

The term I[ containing a self-energy insertion for the other particle of species b has a similar shape 

^ (*.*') = J dt 2 di 2 {g< (t,t 2 ) 9 < (t,t 2 ) [(-*)SM(*2,*2)Sf(t 2j t2)] H)4o(*2.*')^o(*2,0 
to 

+Wa, {tM)g b R fi {tM [gtfi^M s^(* a ,*a)] (-Vgiofat') g^fat') 

+igK (t,t 2 )g b H :o (t,t 2 ) [t^ (t 2> t 2 )Ef (t 2> t 2 )] <£o(*2,0 < (*2,i')} ■ 



(A4) 



The third first-order term in the perturbation expansion is the ladder term. Each of the two vertices can have the 
Keldysh indices + and — . Thus one gets the following four terms 
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l[ 3) (t, t')=ij dhdh [g a>0 {t, h) g h;0 (t, h) Vg{h,h) 9^o(h,t') <? 6 < (i 2 , (A5) 

~9a,0 9b,0 Vab > 9a ,0 9b,0 ~ 9a : 9^0 Kl> < 9a,0 9b,0 + 9a,0 9b,0 ^ab 9afi9b,0 • 

The causal and anticausal Green's functions can be eliminated in favour of retarded and advanced GF's, see Eq. (g), 



I 



(A6) 



-- 9a,o 9b fl Wab 9a \o 9b,o + 9a,o 9b,o Wal R 9a fl 9b fi + 9a \o 9b ft Mai A 9a fl 9b,0 
+9a,o 9m iVab < 9a,o 9b,o + 9a,o d fi W a s b > gio 5b < 

+.9^0 9b fi iVab R 9a fl 9 b fl + 9a fi 9bfi Wab A 9a fi 9b fi + 9a fi 9b fi ^ab 9a fi 9b fi ■ 

There are three classes of terms in the above equation: (i) terms ending with a product 5a,o3&,0' ( n ) terms beginning 
with one retarded function and ending with one advanced function, and (iii) terms beginning with j^gj^. 
The further procedure is presented in detail for the first term in the above equation 



dhdt 2 g* (t, h) gfM, h) iVM^h) g<o(h,t') 9b. (h, 



(A7) 



dhdhigZaihh) < (M~2) V ab S(h-t 2 ) + e(t 1 -h)V^ > (t u t2) 



+e(h - h)v a s b <(h,h)\ g< (h,t') 5 <# 2 , . 

The Heaviside functions allow it to use the semi-group property in certain functions g R and g < , respectively, in the 
following manner 



J dt x dt 2 | 



<fl2o(tJi)fl?o(*.fe)J [VabS(h-h)\ [^oft.O^o&.O 
+ [i9a, {h h) g^oit, h)] [«?£o(*i,*2) ivs>{h,h) 9 R ,o(h,h)] [g< (h,t') g^ Q {h,t') 
+ [ig%{h h) g% (t, h)] [ig R , (h,h) iV a s b <(h,h) s£o(*a>*i)] [g^o(h,t') g^ihJ) 

Renaming the integration variables we arrive at 



(A8) 



dtdt 



+ 



V ab S (t - t) + ig? (i, t) iV a s b > (i, i) g£ (t, i) 



(A9) 



ig% (t, t) iV a s b < (t, t) g R (i, t) 
Thus this term belongs to the anticipated structure 

dUtQX b {t,h\V R ab{lt)Q<b{lt l ) ■ 
(3) 

Making the same analysis for all term or I\ , one gets 

I?\t,t') = [ dMig* (t,hgb R A t ^\ V «bS(t-t) + ig£ (tJ)iV a l > (l 



I- 



(A10) 



(All) 



+ig< (t, t) iV a s b R (t, t) g R fl (t, t) + ig R (t, t) iV a s b < (t, t) g R (i, t) 
+ig< (i, t) iV a s b A (t, t) g R (t, t)}g< (t, t') g< (t, 

- / didtigUt, t) g R (t, t) \g<(t, t) iVf b <(t, t) g<(t, t) 



1G 



-g< (t, t) iVf b > (t, t) g< Q (i, t) {-i)g% (f , t>) g£ (t, ^ 



dtdtg< (t, t) g< (t, t) lV ab S(i- t) + (-i)g£ (i, *) t) g£ (i, t) 

+ H)9a.o(t, t) V a s b R (t, t) ig< (t, t) + (-i)g£ Q {l f) iV a s >(t, t) g£ (i, t) 
H-i)9b,o& *) Vf b A (i, t) ig<(t t i)\(-i)gUt, t') gUl t') . 



A comparison with the structure ( p6| ) gives 4 additional terms for E^ b (S^ b ), and two further terms of a< b . Alltogether 
we get the expressions (^l|) for E^and ([||) for a< b . 



2. Second— order contributions 



According to Eqs. (^6|) and ( |37| ) the second-order terms should lead to diagrams with two self-energy insertions of 
first order with respect to V s (reducible diagrams) as well as to diagrams with one sclf-cncrgy insertion which is of 
second order. We will demonstrate this here for one typical term. The analysis of the ladder term with two rungs (cf. 
Fig. |) leads, among many other terms, to the following contribution (ii, £2,^1,^2 are integration variables) 



ff"o(*. *i) S*o(*> fc) iV^(t x ,t 2 ) 9%{hM) 9 R (t2, h) (A12) 

x iV^hM) 9< (h,t') gUh,t') . 

The procedure to achieve a two-time structure is similar to that in t he fo regoing subsection. According to (jll^) each 



causal function V s consists of 3 terms what leads to 9 terms in Eq. (A12). All contributions containing at least one 



time-diagonal part are easily shown to be reducible. Therefore we concentrate on the others 

I'z = J 9a,S, ti) <$>(*. *a) Wi - t2)V>(t 1 ,t 2 ) + 0(t 2 - ti)V„1(ti,fc)] 9a fi (h,ti) ggafah) 

- h)V> b {hM + 6(t 2 - h)V< b {hM9to(ti,t')g^ {h,t') ■ (A13) 
These terms should fit into the following structure 

j didig* b {t, t) (2) (i, t) g<(t, (Au) 

+ J dtdi x dtidtG* h (t,t) E&<i)(Mi) G* b {h,ti) Xab(i)(tui) Gab&t') , 

where 2^(2) denotes the two-particle self-energy in second order, and S^q) are the first-order quantities identified 

in the foregoing subsection. 

Analyzing the expressions in Eq. ( A13| ), we find that the "mixed" terms (with one V < and one V > ) are reducible. 



The terms containing two functions V^(pr two functions V^) lead to a reducible term for t 2 > i\ (t\ > i 2 ) and to 
an irreducible part for i 2 < i\ (t\ < t 2 ). This is shown in Fig. || in form of diagrams. The reducible terms contain 
two-particle self-energy insertions of first order with respect to V . The second-order terms contributing to £^,( 2 ) 
are given by 

£& (a) (t,0 = / dhdt 2 e(h -t 2 )55(*,*2)v;>(t ) t 2 ) 5 5(i,ti)^ (t2,Ov^(ti,Ofl5( i i.*') 

+[ll'a< — >22'b\. (A15) 

This term can be shown to be a vertex correction to the two-particle vertex. 

The other second-order diagrams in Fig. ^ can be discussed in a similar way. The terms with two single-particle 
self-energy insertions for the same particle are reducible in any case. For the other two types of diagrams, there are 
reducible as well as irreducible parts. This is shown in Fig. |[ The first diagram in each row is a reducible one. The 
second diagram is not reducible and it corresponds to a vertex correction term. The third diagram is not reducible 
as well, but is the first self-energy correction to the diagrams of the two-particle self-energy of first order (cf. Fig. ||). 

Because we started from a ladder equation we do not find all possible second-order terms contributing to the two- 
particle self-energy. Therefore we will restrict ourselves to the self-energy in first order with respect to the dynamically 
screened interaction. 
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3. Analysis for the other correlation functions 



The analysis for the functions g^ b + , g ab ++ , and g ab ++ can be made in the same way as above. It is scetched 
here for <?+ fc h . Consider the first rung diagram. Evaluation on the Keldysh contour gives in analogy to Eq. ( |A6| ) 



r(3) _ R R - V S < > , R < - V S R < > , < R -y-S 'A < > 
1 1 — 9a,o9b,O lV ab9a,()9b,0 + 9a,o9b,O lV ab 9a,0 9b fl + 9a,0 9b,0 1 V ab 9a,0 9b,0 



(A16) 



+9a,0 9^,0 iVab 9a,0 9b,0 + 9a,0 9b,0 W ab 9a,0 9b,0 
+.9a,0 3^0 iVab R 9a fi 9b,0 + 9a,0 9b,0 ^ab A 9a 9b + 9a,0 9b,0 ^ab 9a fl 9b,0 ■ 

Again the two-time structure can be achieved, and the structure is (cf. Eq. |3q) 



9ab 



r R 

y a b 



V ab 



^ab a ab + ^ab + Gtf, 



^ ab 



Vab 



y a b 



(A17) 



with Q^ b + = g<o g b0 - The fourth and the fivth term in ( A16) are contributing to a^ b + the following: 9h > Vf b g^ + 
9a0^ab9b > o- Together with the respective contributions from the diagrams involving single-particle self-energies, 



a ab i s gi ven by Eq. ( |5l| ) then 



APPENDIX B: KRAMERS-KRONIG RELATION FOR THE TWO-PARTICLE SELF-ENERGY 



Hermitean and antihermitean parts of the self-energy can be constructed in the usual way by 
E™ (121'2' ; uT) = l - [Ef b (121'2' 5 loT) + S« (1'2'12, ujT)*] , 
S^(121'2', uT) = l - [E« (121'2' 5 loT) - E* (1'2'12, wT)*] . 



Due to Ef & (121'2',wT) = E^(l'2'12, loT)* we have 

E« 6 (121'2', loT) = 1 [s* (12l'2', W T) + E^(121'2', wT) 
E^(121'2', loT) = 1 [e* (121'2', W T) - E^(12l'2', wT) 

Because of the causality there exist the following Kramers-Kronig relations (see, e.g. 

du>' E£(u/T)-£& (T) 



E&«T) - E™(T) = *P /" ^ 

J-oo 71" 



and 



S^T)-E^°(T) = iP /"^ 

J-oo 7T 



du;' E»(u,'T)-E*°(T) 



where P denotes the principle value. The static parts, 

1 



E«°(T) = Ef F + Ef F + -(7V ah - 1)K 6 + ^K 6 (AU - 1) 
s^°(T) = -W b - i)y a6 - \v ab {N ab - 1) : 



(Bl) 



(B2) 



(B3) 
(B4) 

(B5) 
(B6) 



have to be substracted in order to ensure a proper behaviour in the frequency integrals. 

In the approximation which is used in this paper both contributions to the two-particle self-energy, A ao and V° b , 
can be split into hermitean and antihermitean parts. For A a & there follows 
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A a H b (pip 2 , p'iP' 2 , n) - (27r) b J pi , pi 5 P2 , p , 



Ir ,_N l±/a,p 1+ q + ns(^l) 



d 3 q 



7T 



" — ^1 — £a,pi+q — Cfc, 



■P2 



and 



(2tt)3 

+ (a < ► b . I — 2i J- . 

d 3 r 



(B7) 



^3 Im e K (q, ft - e a:Pl+q - e b , p 
xsKa(q) 1 ± / ajPl+q + n B (si - e a , Pl+q - e b . P2 ) + (a < — — >2) 



(B8) 



Because the matrix elements of A are diagonal for homogeneous systems, the hermitcan part is a real quantity whereas 
the antihermitean part is imaginary. 

Hermitean and antihermitean parts of V^f are given by 



Kf W (PiP2,PiP2,") = (2^) 3 <5 Pl+P2 , pi+P2 V o6 ( Pl - pi) 



(B9) 



7T 



Ime R (pi-p'i,wi) 



1 + n B (vi) ± \fa, Pl ± \ fb,y' 2 



±1 



+i Im e R (pi - pi , - e 0)Pl - e b<p > ) — / QiPl - / 6)P 



+ (a< — ► 1< — >2) 



Kf ^(pip 2 ,pip 2 ,0) = (27r) 3 ( 5 Pl+P2 , p;+P2 y ob ( Pl - pi) 



duj\ 
n 



- lme R (pi - pi, wi) 



-z line 



+(a 



Pi -pi,^-£a, Pl -e 6 , pi J(^l + nB^O-e a , Pl - £b,p 2 J ± g/a.Pi ± q-^.p- 



6, 1 



(BIO) 



Here, the matrix elements are complex in general. Diagonal matrix elements in an arbitrary representation, (n|V a c 
and (n.|V^jf \n), are of course real and imaginary, resp ecti vely. 

It is easy to see that the Kramers-Kronig relations, (B3), are fulfilled. 



off Hi 
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APPENDIX C: CONNECTION OF Gab TO THE MATSUBARA GREEN'S FUNCTION 



It is the purpose of this appendix to scetch for the case of thermodynamic equilibrium the connection to the 
Matsubara technique for imaginary-time Green's function. Especially, it will be shown that the Fourier transform of 
the advanced two-particle Green's function, which was defined in Eq. is closely related to the analytical continuation 
of the two-frequency Matsubara Green's function. 

Imaginary time Green's functions dependend on three times were considered by Rajagopal et al. in connection 

with an Bethe-Salpeter equation for an electron-positron system. In this appendix the same notation is used. A 
function in which the times of the creation operators are equal is given by 



G^ihhtz) = -{T w (^(ti)</>(ta)B{t a ))) 
i 



(CI) 



with 



B(t 3 ) = ^(t 3 )^(t 3 



(C2) 



This function depending on three imaginary times is just the function for which the Bethe-Salpeter equation is 
"closed" . The term "closed" means here that both sides of the equation contain the same type of functions. 

For the Fourier coefficient in Matsubara representation, Rajagopal and Cohen |^2j found the following important 
double dispersion relation 



with 



U>1 — U>1 



/ (1) (^1^2) 



U>\ + UJv — UJ\ — Ll>2 



1 



P 

UJi> — U>2 



(C3) 



/ (2 W 2 ) 



d(h - fa) d(i a - is) e --i(*i-*3)-Wt 2 -ta) _/ ^ (tl) [ B (t 3 ),tf>(t 2 )]- 

i \ L 



(C4) 



and 



fW{w lU2 ) = r d(h ~ t 3 ) d(t 2 - t 3 ) e -^{t,-t 3 )-^{u- t3 ) l^tt), cf>(t 2 )}+ , B(t 3 )] _ ) 



(C5) 



The Matsubara frequencies are cj" = (2v + (3 + ijS a \ 

Assume that the function & ep ' (lu^ujP) has been determined. Then the question arises how to extract from this 
quantity the physical information on the two-particle problem. The usual way is to carry out the sum over one 
frequency, 



(C6) 



and performing afterwards the analytical continuation. The summation, however, leads to the Matsubara coefficients 
of the causal Green's function depending on two imaginary times, cf. (|l9|). 

We want here to proceed on an other way. The analytical continuation is done for the two-frequency quantity. This 
function G^ ep '(zi, z 2 ) is considered now for zi = Q e — irj and z 2 = fl p — it] with 77 — * 0. Making the transformation 



f2 e , fi p — > f^ep, we have 



1 



U2) 



ep 



Q p — us 1 — ii] Q ep — uj\ — UJ2 — ifj 



(C7) 



f2 p — ui2 — if] 



Now an integration is perform with respect to £l p what leads to (aditionally the variables lj\ , 0J2 are transformed into 



20 



oo oo 



/ (2) (CJ12 L0 2 ) / (1) (^12^2) 
f2 ep - fip - + UJ2 - ir) &ep — U\2 - IT) 



(C8) 



1 



fip — UJ2 — ir) 

00 



(27r) Slep - W12 - ir] J (27r) 



/ (2) (W12 w 2 ) + (wi 2 w 2 ) 



(C9) 



The integration over w 2 can be performed and gives the Fourier transform of two-time correlation functions 

J ^f(Lo 12 LU2) = J°° d(h - i 3 ) d(t 2 - h) e-*-i»(*i-«s)-^(t a -tx) i2j tg) (CIO) 

d(i 1 -t 3 )e-^ 12( * 1 - t3) /(ti,ti,t 3 ) 



The equal-time commutator [ip(ti), 0(ii)]+ vanishes and therefore also the contribution with /Win (C9). The final 
result is 

00 



7 (27r) 12 ep - W12 - j?/ 



with 



f^M =ld(t- t') e ---(*"*') i( [^(t), [^(f')^(f'), 0(i)]_ ] + ) . 
This result has to be compared with the Fourier transform of G^ b defined in (|2 



du)\ 



Gabi^ab) — / JZ \ ^ 7 

J (2tt) \l ep - u 12 - 1 



hab(uJl 2 ) ■ 



with 



h ab (u 12 ) = fd(t- e -*"«<*-*') [* (n, t), [*„(r 2 , t), *J(r' 2 , O^ft, 
This is, except a minus sign, the same expression. 
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FIG. 1. Time ordering on the Keldysh double time contour for the function g£ b ' • The ordering is causal on the upper 
branch and anticausal on the lower branch. 
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FIG. 2. Diagrammatic expansion of the dynamically screened ladder equation up to second order. 




FIG. 3. Diagramatic structure of the two-particle self-energy in first order in V s . 
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increasing time i 




■ i i 

■ i i 




FIG. 4. Evaluation of the ladder terms with two rungs. The first, the third, the fourth, and the fifth term are reducible, 
i.e., there are two successive self-energy insertions of first order. The second and the last term, however, are not reducible: they 
are two-particle self-energies of second order in V 5 . 




FIG. 5. Here two other types of second-order terms are shown: in the first row two single-particle self-energies are combined, 
in the second row one single-particle self-energy is combined with one interaction between the particles. In dependence how 
the times do overlap, there are three different kinds of diagrams. 
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